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1. Introduction 

This paper develops a tlieoretical approacli allowing one to consider PCI in the Wannier 
regime. The interest in this problem has been initiated by recent experiment |I| followed 
by theoretical work |0] in which the electron impact in the vicinity of the atomic inner 
shell ionization threshold 

e + A ^ e^ + e^ + A+* (1.1) 

was studied. In this reaction the inelastically scattered electron e^ and the knocked out 
electron e^ have low energy propagating in the Coulomb field of the positive ion A'^*. 
These low-energy electrons are known as the Wannier pair. The Coulomb interaction 
of Wannier electrons between themselves as well as with the field of the ion results in 
the Wannier power law for the cross-section 

a^e^" (1.2) 

where e is the energy excess over the ionization threshold and the exponent /x depends 
non-trivially on the ion charge Z 

VlOOZ - 9^ 
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A distinctive property of reaction (|1.1|) is the finite life-time of the inner-shell vacancy 
in the ion A+*. Its decay 



e„ + e/3 + A+* ^ e„ + e/3 + e^ + A2+ (1.4) 

creates the fast Auger electron e^ and doubly-charged ion A2+. The fact that the pair 
e^, e/3 has low energy guarantees that at the moment of the decay it is situated in 
some close vicinity of the ion. Therefore the pair should strongly interact with the 
products of the Auger decay, i.e. the doubly charged ion and Auger electron. An 
interaction of this type is known in literature as the Post- Collision Interaction (PCI). 
What distinguishes reaction ( p..l|) , (|1.4|) is the fact that the Wannier pair takes part in 
PCI. Previously PCI was extensively studied in the near threshold region for the case 
when a single low-energy electron is localized in the vicinity of the atom at the moment 
of the autoionization decay. 

Thus reaction (|1.1|) , ([T7^ ) combines the two problems of the final state interaction, 
the Wannier problem and PCI. Previously they both had been extensively studied 
experimentally as well as theoretically but in those situations when only one of them, 
either the Wannier problem or PCI had manifested itself in the final state. 

The main purpose of this paper is development of the theory describing PCI in 
reaction ( p. . 1| ) , ( p^ ) as well as in similar situations in electron or atomic collisions or in 
photoionization. One of the interesting qualitative questions arising from the reaction 
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( p..l[ ),( p!^ is a possible influence of PCI on the Wannier power law ( |1.2| ). In the field of 
the singly charged ion the Wannier exponent ( |1.3|) is ^ 1.127, while in the field of the 
doubly charged ion it is /i2 — 1.056. If the energy of the Auger electron e^ is high enough 
then it does not interact strongly with the pair eQ,,e/3 and the decay could be considered 
as a switch of the ion charge from the intermediate-state value Z = 1 to the final-state 
value Z = 2. The question is how this charge variation affects the exponent. We will 
demonstrate that for sufficiently low above-threshold energy the exponent is equal to 
H2, while for higher energies it is fii. This fact agrees with the following simple physical 
picture. When above-threshold energy is very low then the Wannier pair is located so 
close to the ion at the moment of the decay that its separation is negligible compared 
with an important in the problem radius of the Coulomb zone. Therefore one can neglect 
this separation, neglecting as well the very fact of existence of the intermediate-state 
resonance. Thus in this case the pair is influenced by the field of the final-state ion 
with the charge Z = 2. In contrast for high above-threshold energy the Wannier pair is 
located well outside the radius of the Coulomb zone at the moment of the Auger decay 
and therefore is not strongly influenced by the field of the doubly-charged ion. The 
main characteristics of the pair in this case are governed by the ion in the intermediate 
state with Z = 1. 

The paper is organized as follows. Section |^ gives a short introduction to the 
Wannier problem which was theoretically studied in a number of papers. The variety 
of approaches developed and the comprehensive list of references can be found in |^ - 
0. Section ^ presents a short review of essential ideas of PCI. The literature on the 
theory of PCI can be found in 0. These introductory notes permit presenting basic 
ideas of the two problems from a particular point of view which provides the easiest way 
for further development in Sections § and ^ where a way is found to describe PCI when 
Wannier electrons are involved. 

Generally speaking the problem considered looks very complicated being dependent 
on the Coulomb interaction of three electrons and the ion. There are, however, important 
simplifications. Firstly, we assume that the Wannier pair possesses low energy, while the 
Auger electron is so fast that one can neglect its interaction with the pair. This allows 
us to reformulate the problem in terms of the Wannier pair only. Still it is a dynamical 
three-body problem. Another important simplification comes from the cornerstone of 
the Wannier problem 0. In order to escape from the ionic field Wannier electrons 
should move along a particular classical trajectory. This motion can be considered as 
a collective mode describing essential properties of the system. One can easily verify, 
see Section |^, that this mode can be described in terms of an effective single-particle 
Coulomb problem. Of course there are important corrections arising from the three- 
body nature of the problem and making it so specific and interesting. The point is 
that knowing the reliable initial approximation one can find these corrections in a clear 
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analytical form. 

These arguments show that one should expect the problem considered to have a 
clear solution. Finding and describing it in simple physical terms is the major goal of 
this paper. The main result obtained is the formula which presents the autoionization 
profile -Rpci distorted by PCI in the form 

i?pei = K^Rc (1.5) 

The factor R^. here arises from the most simple approximation for the Wannier problem 
described by the above mentioned effective single-particle Coulomb problem. Due to 
this reason we will call it the Coulomb factor. This factor depends on the overlapping 
integral of the Coulomb wave functions. It is interesting that a similar factor describes 
PCI in the traditional situation when there is only one low-energy electron. The only 
distinction is that in the later case there is a true single-particle Coulomb problem, while 
for the Wannier pair the Coulomb problem arises as an effective approximation. Still 
the similarity is so close that one can consider the Coulomb factor Rc as a well known 
quantity which was well studied previously both experimentally and theoretically. In 
contrast the factor K„ takes into account more subtle and interesting properties of the 
process which are essential for PCI with the Wannier pair and have no analogy in PCI 
with one low-energy electron. Due to this reason we will call it the Wannier factor. In 
Section § an explicit analytical formula for this factor is derived. The Wannier factor 
found depends on the normalization coefficient of the Wannier electrons wave function. 



It is interesting that the Wannier power law ( |1.2| ) arises from the properties of the 
same coefficient considered in the region of small separations of the pair from the atom. 
In contrast, when PCI is studied this coefficient manifests itself in a wide range of 
distances. This property of the profile ( p..5| ) opens a possibility for experimental studies 
of the normalization coefficient of the Wannier pair wave function in a wide area of 
distances. 

The results obtained in Sections are illustrated in Section || for the electron 
impact near-threshold K-shell excitation of Ne. Section ^ presents a short description 
of the major results. 



2. Wannier problem 

Consider two low-energy electrons in the field of an ion with the charge Z. The Coulomb 
repulsion between the electrons makes opposite directions the most favorable for them 
0]. To simplify the following consideration let us use this fact developing the following 
model. Let us assume that the electrons are moving in opposite directions neglecting 
fluctuations in their angular distribution. Moreover, let us neglect as well the orbital 
moment. Thus in this model angular variables are neglected which makes it essentially 
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one-dimensional. This model was applied earlier for description of the Wannier power 
law It is important that the most interesting part of the Wannier problem, its 

instability against a non-symmetrical distribution of energy between the two electrons, 
is taken into account in full. Moreover, we will argue in Section ^ that our final results 
do not depend on model assumptions. One of the main advantages of the model is 
simplification of notation, which inevitably is to be sufficiently complicated. 

Our strategy in this Section is to first consider the classical approximation and then 
develop the quantum picture using the semiclassical approach which gives a correct 
description of the process for low above-threshold energy. The semiclassical approach 
to the Wannier problem was developed in detail in to describe the near-threshold 
behaviour of the reaction ( |1.1|) . In this case one is mainly interested in the behaviour of 
the wave function of the Wannier pair for small separations. We will see that to describe 
PCI it is necessary to know the wave function for any separation, not necessarily in the 
vicinity of the atom. The convenient representation for this wave function valid for an 
arbitrary separation is given below in equation ( |2.62|) which is the main result of this 
Section. 

2.1. Classical trajectories 

Consider classically the problem of two electrons e^e^j in the framework of the above- 
formulated one-dimensional model in which electrons are supposed to move in opposite 
directions with zero orbital moments. Such a system is characterized by two radii 
f'ayfis > 0. The Lagrangian in these coordinates is 

L{r^, rp,r^,rp) = ^-f + '^ + - + - (2.1) 

2 2 r^ rfs r^ + rfs 

where Z is the charge of the ion. It is convenient to introduce hyper-spherical coordinates 

(2.2) 



r„ = r cos(x + f ) 
rp = rcos(x - f) 



in which r > describes the overall distance of the pair from the ion, and —it /A < 
X < 7r/4 shows how strongly the pair is deviated from the symmetrical configuration 
r^, = r^, in which x = 0. In these coordinates the Lagrangian is 

• 2 2-2 

Hr, r, x) = y + ^ - f/(r, x)- (2.3) 
The potential energy 

U{r,x) = -^ (2.4) 
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can be looked at as a Coulomb field produced by a x-dependent charge 

, , (2Z - l)cos2y + 2Z , , 

^ix) = —7T^ > 0- (2-5) 

v2cosxcos2x 

The function w{x) has a minimum for the symmetrical configuration x = 0- It manifests 
itself as a maximum at x = of the potential f/(r, x) considered as a function of x- This 
property of the potential is known as the Wannier ridge. The ridge makes the classical 
system unstable. The instability results in the most important property of the Wannier 
problem. In order to escape from the Coulomb field of the ion the two electrons should 
propagate in the vicinity of the symmetrical configuration. Otherwise one of them would 
be trapped by the ion field, i.e. the pair fall off the ridge. That is why the motion in 
the vicinity of the top of the ridge, where the pair is stable plays such an important role 
and is considered in detail below. 

In the vicinity of x = the function w{x) can be expanded in the Taylor series 



w" 



w{x) ~ w + — X , (2.6) 



where the coefficients are 

AZ -I , , 

» = (2.7) 

„ 12Z - 1 , , 

(2.8) 

Equation (|2.3| ) results in equations of motion 

re-^, (2.9) 



■2 ^(X) 

A. 

d , r, 1 dw 



i^x) = ^X)- (2.10) 
The trajectory along the top of the Wannier ridge satisfies 

Xo = 0, (2.11) 

fo = - ^. (2.12) 

The second equation coincides with the equation describing the motion of a single 
particle in the field of the Coulomb charge w. We will refer to this motion, which 
plays an important role below, as an effective single-particle Coulomb problem. The 
first integral of (2.12) 

f-- = ^ (2.13) 

describes the energy conservation law. Using it one finds the ridge-top trajectory 

1/2 



+ 



w 



2e 



2ero(t) 



w 



1 + 



w 



1/2 



For small deviations from the ridge-top trajectory one obtains the linear equation for 
X = 5x using expansion (|2.6|) 



(2.15) 



To simplify notation the ridge-top trajectory is denoted here by r = ro(t) defined in 
(2.14). For further applications we need to consider solutions of (|2.15|) in detail. It is 
convenient to change the argument t — x in such a way that 
d w d 



dt 

This is achieved for 
1 

^=2 



2^/2edx' 



1 + ^ 

er 



1/2 



+ 1 



The new variable x = x{er) satisfies 

,1/2 



X 



1 / w 

2 ler 



1 + — 



er < 1 
er > 1 

Using (|2.16| ) one can rewrite ( p.l5|) as 
d'^xix) A 



dx'^ x{x — 1 
where A is the coefficient 



VIZ 



X = 2— = 2- 

w AZ -1 



Scaling the function 

X{x) y{x) 



X{x)/x 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



we find that the motion in the vicinity of the top of the ridge is described by the following 
linear differential equation 



x(l - x)y"{x) + 2(1 - x)y'{x) + Xy{x) = 0, 



(2.21) 



in which one recognizes the known [|TT| hypergeometric equation whose parameters 



a = u + 1, P = -u, 7 
are restricted by the condition on u 
u(u + l) = A. 



(2.22) 
(2.23) 



Solving this equation and using ( |2.19| ) one finds an explicit form for this parameter 
"/100Z-9\^/2 



1 

z/ = - 
2 



\ AZ -1 



1 



(2.24) 



which governs the Wannier problem, compare ( p.24|) with ( |1.3|) . 

It is convenient to choose a set of basic solutions yi{x),y2{x) of ( p.21| ) as 

yi = ix- ly-'^-'^F (7 - a, 7 - /5, 7 + 1 - a - /?, 1 - x) , 
?/2 = x~"F (a, a + 1 — 7, a + 1 — /?, 1/x) , 



where -F(a, /3, 7, x) is the hypergeometric function. Then an arbitrary solution of (|2.1^ 
can be presented linear combination 



x{x) = aif{x) +a2g{x) 
of functions f{x) = xyi{x), g{x) = xy2{x) having the following form 
f{x) = x{x - 1)F {-V + 1, z/ + 2, 2, 1 - x) 



9{x) 



c 



Z/+ 1 
2v + l 



x-T (u,u + l,2u + 2,-] . 



X 



(2.25) 

(2.26) 

(2.27) 



To simplify the subsequent formulae a normalization coefficient ( in (2.27) is chosen to 
be 



C 



(2.28) 



r(2z/ + 1) ■ 

Here T{x) is the usual Euler gamma function. We will need to know the asymptotic 
conditions for the functions f{x),g{x), which are easily found from (2. 26), (2. 27) 



1 



C{^+i) 

X — 1 



X 
X 



00 

1, 



(2.29) 



g(x)^\ (2.30) 
[ 1 X — > 1. 

At the end points x = 1 and x = 00 one of the functions, either /(x) or g{x), is singular 
while the other one is regular making them convenient for applications. Notice that the 
derivative g'{x) remains singular at x — 1 

w" 

g'{x) 2— [ln(x - 1) + d] + 0[{x - 1) ln(x - 1)]. (2.31) 
Here c? is a constant 



d = 2 [^{u + 1) + C] , 



(2.32) 
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where iplx) = r'(a;)/r(a;) and C = —ipil) — 0.577 is the Euler constant. One more 
necessary parameter in the problem is the Wronskian of the functions f{x),g{x) which 
can be found from asymptotic conditions (2. 29), (2. 30) 

f'{x)gix) - f{x)g\x) = 1. (2.33) 

The given consideration describes in detail the main trajectory along the top of the 
Wannier ridge (2.12) as well as small deviations from it given in (|2.25| ),(2.26,(2.27) in 
which the variable x = x{er) is defined in ( |2.16| ),(2.14). 



2.2. Quantum description 

Essential to the problem distances are in the Coulomb zone 

r ~ i = rc. (2.34) 

In this region the Coulomb force is Fc ~ ~ e^. The variation of the Coulomb 
energy on the wavelength {l/^/e)Fc ~ e^^"^ proves to be smaller than the kinetic energy 
~ e provided the energy is low 

e<l. (2.35) 

It is supposed in this paper that the later condition is valid. The given estimations show 
that in this case the semiclassical approximation turns out to be correct permitting one 
to write the wave function in the well known form 

M(l)=N{q,P)exp[zS{q,P)], (2.36) 

where S{q, P) is the classical action. The wave function depends on a set of coordinates 
q and quantum numbers P, the later ones play a role of parameters. It is convenient to 
choose P as momenta conjugate to coordinates q. To be specific consider the convergent 
wave function for the final state ip = ip^~^ which is the most interesting in the Wannier 
problem. This wave function may be looked at as a transition amplitude from the state 
in which coordinates q of the Wannier pair are in the vicinity of the atom to the final 
state in which momenta P describe the motion at an infinite separation. In classical 
terms the transition should be described by a trajectory with initial coordinates q and 
final momenta P. These two sets of parameters q and P fix the trajectory. The classical 
action depends on both these sets of variables S = S{q, P) . It is well known that the 
non-exponential factor N{q, P) in the wave function depends on the derivatives of the 
action |T^ 



N{q,P) 



det 



dqidPj 



1/2 



(2.37) 
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Let us choose a set of coordinates g as g = (r, rj), where rj is the momentum conjugate 



to the variable x introduced in (2.2) 

V=jr = 2.38 
OX 

A convenience of this choice becomes apparent soon. The corresponding momenta are 
P = {p,x)- (More accurately the momentum conjugate to rj should be written as — 
The minus sign is dropped to simplify notation.) Here p is a momentum conjugate to r 

dL 

P=^ = r. (2.39) 

To avoid confusion the bar symbol is used to mark the variables belonging to the final 
point of the trajectory for infinite separations. In this notation the energy of the pair is 
written as e = p^/2, and the action looks like S = S{r, ri;p, x)- Calculating the action 
we should consider the region in the vicinity of the Wannier ridge where both x and rj 
are small. Therefore we can expand the action in powers of these variables 

S{r,v;P,x) = So{r,p) + f:{r,ri;p,x), (2.40) 

where 

S'o(r, p) = S{r, 0; p, 0) = £ ^p"^ + 2^dr' - pf (2.41) 

describes a one-dimensional motion along the ridge-top trajectory governed by the 
effective Coulomb charge w. We will recognize later, see ( |2.55| ), that a slightly different 
choice of the effective charge proves to be more convenient. The quantities which will be 
recognized be affected by this redefinition are marked by the hat symbol. The quantity 
S(r, r7;p, x) depends on the derivatives of S. The first derivatives for the ridge-top 
trajectory are obviously zero. Therefore, the expansion starts from the second powers 
of x^rj. It is shown below, see discussion after ( p. 621 ), that it is sufficient to restrict our 
consideration to these lowest order terms writing 

— 2 2 — 

S(r, m P, X) = ^A + xvB + ^C. (2.42) 

The factors 1/p and p are included here to simplify the subsequent formulae. The 
coefficients of this expansion are functions of r and p. In order to find them let us 
remember that x and 1] are the conjugate variables and therefore 

li - ^ - (2 A3) 

P dx^ \9xJr,' 




pC 



drf 
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(2.45) 



The derivatives in ( p . 43| ) , (P . 44|) , ( |2 . 45|) describe the variation of the coordinate and 
momentum under variation of the trajectory. This variation in turn is described 
by variation of the coefficients ai,a2 in (|2.25|) which governs the coordinate x{^) = 
aif{x) + a2g{x) as well as the momentum rj conjugate to x 



w w 

= - — (ai/'(x) +a2g'{x)). 



(2.46) 



Using (^^),{^^ one finds from (^) 

{dai)f' + {da2)g' 



2 



{dai)f + {da2)g 



J ri 



where according to ( |2.46| ) the condition rj = 
dai, da2 

{dai)f' + {da2)g' = 0. 
Combining the last two equations we find 

f'g' - f'-g' 

2 f'-g - fg' ■ 

The other derivatives are calculated similarly 

f'g - fg' 1 



const restricts variations of coefficients 



B 



C 



f'g-fg' f'g-fg''' 



2 f-g - fg 



(2.47) 



(2.48) 



(2.49) 



w f'g - fg' 

Evaluating the last identity in ( p.48| ) the Wronskian ( p.33|) was used. All three 
coefficients turn out to be functions of the single variable x, A = A{x), B = B{x), C = 
C{x). This property of the coefficients was anticipated when normalization factors 1/p 
and p were introduced in their definitions in (|2.42| ). 

We can now use the fact that for the final state f = oo, and therefore according to 
1. From (2.29), (2.30) we find / = 0, 

w f g'{x) 
2 



/' = 1, (? = 1, resulting in 



A{x) 



fix) 



g 



B(x) 



C{x) 



f'ixy 

2 fix) 
w f'{x) 



(2.50) 
(2.51) 
(2.52) 
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All these quantities are well defined functions of x, except for the constant g' which 



according to ( 2.31 ) is logarithmically divergent g' 
divergence manifests itself in the action term 



ln(x — 1), 



X 



oo. 



This 



wx 



4p 



-9 [x) 



2p 



In 



2p^r 



w 



d 



(2.53) 



In order to deal with this problem let us remember that the Coulomb action 5*0 ( 2.41| ) 



reveals the logarithmic divergence as well, 5*0 ~ {w / p) \n{2p f / w) , f oo. This 
later divergence proves harmless. Really, it takes place for an effective one- dimensional 
Coulomb problem described by the only coordinate r. For this simple situation 
the Coulomb logarithmic divergence is of course well studied and cured by textbook 
methods, see [O]. Having this in mind let us "renormalize" 5*0 and S, 5*0 Sq, E — E 



Sir,r];p,x) = So{r,p) + T.{r,r];p,x) (2.54) 
in a way that makes the divergence (|2.53| ) apparently harmless as well. The action 



5*0 (r,p) in (|2.54| ) is defined as 



So{r,p) = Jp^ + 2^^dr' - pf 



(2.55) 



It describes the one-dimensional motion in the field of the Coulomb charge w(x) ( |2.5| ) 
which explicitly depends on x^, see ( p.6|) . This property distinguishes it from the "non- 



renormalized" quantity 5*0 (|2.41| ) which does not depend on x- The renormalized E 



w"x' 



dr' 



r r'^p"^ + 2wlr' 



(2.56) 



is defined in such a way that it compensates for the dependence of 5*0 on x^ . Using 
these new definitions we find 



— 2 2 — 

S(r, Tj] p, x) = —A{x) + xvB{x) + ^C(x), 
2p 2 



(2.57) 



where B{x),C{x) remain equal to definitions given in (2. 51), (2. 52). In contrast A{x) 
differs from the "non-renormalized" value A{x) (2.50), namely 

w g'{x) 



A{x) 



w"h{x). 



2 fix) 

Here h{x) is a new function whose definition follows from ( p.53| ),( p.56D 

dr' 



(2.58) 



h{x) = —p J 



r r'^p"^ + 2w/r' 



w 
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Calculating the integral and taking the limit r ^ oo we find 

r 



h{x) = In 



n w _ / „ 2w 

p -\ \-p\ P H 

w \ r V r 



— (i = In 



X 



X — 1 



(2.59) 



where d was defined in ( p.32|) . 

The non-exponential factor N{q, P) in the wave function ( |2.36| ) is easily calculated 
because with the necessary accuracy it is 



/ \ 1/2 

N{r,r];p,x) ^ No{r,p) ( 1 = No{r,p)B{x). 



(2.60) 



Here the first factor NQ{r,p) is related to the effective one- dimensional Coulomb problem 



No{r,p) 



1/2 



(2.61) 



^ drdp J 

The factor B{x) was calculated previously in (2.51). 

Combining the above calculated quantities we can present the wave function (|2.36| ) 
in the following final form 



(2.62) 



We need no more to distinguish by a bar symbol the variables at the final point supposing 
that p = p, X = X- Evaluating ( p.62[ ) it was taken into account that No{r) exp{iSo) 
gives the semiclassical representation for the Coulomb wave function (j)^~\r,p; w{x)) 



(p^ \r,p;w{x)) = No{r) exp{iSo) 



(2.63) 



This single-particle wave function gives an effective description of propagation of the 
Wannier pair near the top of the Wannier ridge. The propagation is characterized 
by the momentum p and the effective Coulomb charge w{x)- It is often more 
convenient to work directly with the Coulomb wave function rather than with its 
semiclassical representation (|2.63| ), though one can keep open a possibihty to return 
to the semiclassical description when necessary. 

The quantity I]{r,ri;p,x) in (|2.62|) is defined by expansion ( |2.57|) , the coefficients 
A{x), B{x), C{x) of which are given in ( |2.58| ),(2.51) and (2.52). From ( p.62| ) we conclude 
that typical values for both x ^^nd t] are restricted by the available above-threshold 
energy. For x we find ~ P- Similarly, for rj we find that f]"^ ~ prjQ. Here rjQ is a 
typical classical value of the momentum t] which according to (|2.46| ) is estimated as 
Tjo = const /p. For low energy, see ( |2.35[ ), both x and r] are small, ^ Vo- 
This fact justifies expansion ( 2.57| ) of the action in powers of x and rj. 

Along with the wave function ( p.62| ) one can consider another solution for the 
Wannier pair ~ \/B exp{—i'E)(l)^~^\ We use conventional notation which shows explicitly 
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the asymptotic behaviour of the Coulomb wave functions 

(2.64) 



VP 



±i (pr H In 2pr 

\ P J 



The normahzation coefficient 1 / ^/p is convenient for the following calculations, see (|5.2| ) 
in Section |^. At this stage an advantage of representation in which x is a parameter 
and ?7 is an argument of wave functions becomes obvious. The potential f/(r, y) ( |2.4| ) 
approaches the value —w[x)lr for r — * oo resulting in the effective Coulomb charge w[x) 
which in this representation does not depend on arguments of the wave function and 
may be considered as a given constant. We have used this fact above to satisfy boundary 
conditions for large separation. In other representations, for example in representation 
in which x is an argument of a wave function and is a parameter it is more laborious 
to satisfy these boundary conditions. 

For effective use of the wave function (|2.62|) we will need to know the asymptotic 
form of the coefficients A{x\B{x\C{x) in ( p.57| ). It can be found from asymptotic 
conditions (2. 29), (2. 30). Inside the Coulomb zone er -C 1, we find 



A[x) ~ ~w"d, (2.65) 
B{x) ^ Cx-'' = C , (2.66) 



w 

1 2 1 

X — 

V + Iw V + 1 \wp'^r 

Outside the Coulomb zone er ^ 1 we find 



1/2 



C{x) ^ -T^-x = ( rz:^- ) . (2.67) 



A(x) ~ > 0, (2.68) 

B{x) ~ 1, (2.69) 
C{x) ~ - ^ ^ 0. (2.70) 

Using (2. 68), (2. 69), (2. 70) we verify that the phase in the definition of the wave function 
( p.62| ) is chosen in such a way that for r — > cxo 



^P,x('^'^) ^ ^exp 
y/P 



■I , ^ 

t \rix — pr m 2pr 
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(2.71) 



The discussion above shows that the main features of the Wannier problem are 
described by the effective single-particle Coulomb problem with the effective Coulomb 
charge w(x)- An essential property of the single-particle problem is a phase shift. In 
order to find it let us present it in the following form 



5{x) = Ux) + ^5{x)- 



(2.72) 
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Here the first term is the usual Coulomb phase shift which arises due to the effective 
motion in the Coulomb field of the charge w{x) 

4(x) = argr(^l-^^^. (2.73) 

An additional phase shift A6{x) is produced by the term Il{r,7];p,x) in the exponent 
in (|2.62| ). Remember that this term is part of the classical action. The variation of the 



action when r varies from r = oo to r = in the semiclassical approximation is identical 
to the phase shift. Therefore 

w"d 

A5{x) = - [S(r = oo, r/ = 0; p, x) - S(r = 0, = 0; p, x)] = -^x'- (2-74) 

Evaluating the last equahty (2. 65), (2. 68) were used. The minus sign in front of the 
square brackets in ( p.74|) agrees with definition (|2.64| ) of the function 0(~)(r; w(x)). 



From the asymptotic condition ( p.71j ) we conclude that the found phases (|2.73|) , (|2.74| ) 



are the only ones which contribute to the phase shift (|2.72|) . 



The found representation for the wave function of the Wannier problem ( p. 62 ) has a 



clear physical meaning. The Wannier problem takes place in the vicinity of the Wannier 
ridge. This allows to present essential properties of the two-electron problem by the 
effective single-particle Coulomb problem. The Coulomb wave function w(x)) of 

this effective Coulomb problem gives the simplest approximation. The nontrivial three- 
body nature of the Wannier problem is taken care of by the quantities S(r, i]; p, x), B{x) 
which accounts for small fluctuations in the vicinity of the ridge-top trajectory. These 
fluctuations result in two interesting properties. Firstly, the state of the Wannier pair 
on the top of the ridge is unstable in respect to sliding down into configuration when one 
electron is close to the atom while another one is far apart. This reduces the probability 
for the pair to survive on the top of the ridge. This effect is taken into account by the 
normalization coefficient B{x) < 1. Secondly, there appears an additional important 
contribution to the phase shift A6{x) in the effective single-particle problem. 

If the wave function ( |2.62D is used in the traditional Wannier problem to calculate 
the cross-section of a process with two Wannier electrons in the final state then the 
normalization coefficient B{x) should be estimated with the help of (2.66) in the vicinity 
of the atom. This gives \ip\ ~ e'^^^, and therefore the threshold behaviour of the cross- 
section is estimated as a ~ l^/'p ~ e'^^'^ = reproducing the Wannier threshold law 

r2D. 
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3. PCI in the near-threshold region 

The resonant nature of reaction ( |1 . 1| ) , ( P^ ) permits one to present the matrix element 
describing it in the following general form 

(_mmK^. (3.1) 

The wave function I'^i) describes the initial state e + A, |\E'j) describes the set of 
intermediate states eo, + e^ + A"*"* with the fixed resonant state of the ion A"*"*, and 
l^f) describes the final state Cq, + C/j + e^ + A^^ of the reaction described by ( |1.1| ),(|1.4| ) . 



Summation in (|3.1| ) takes into account all possible states for Wannier electrons CcC/j in 
the intermediate state. The potential Vcoi describes the part of the Hamiltonian which is 
responsible for the collision, and Va is the potential which causes the Auger decay. The 
energy denominator depends on the above-threshold energy as well as on the energies 
of the Wannier electrons in the intermediate state and the total width F of the 
autoionizing atomic state. We are interested in the low near-threshold energies for the 
Wannier pair (p.35|) . In contrast the autoionized electron is sufficiently fast and we can 
assume that its energy is so high that 

e < K e^. (3.2) 

The width F is another small parameter of the problem. We will assume that the width 
is much lower than the above threshold energy 

F < e. (3.3) 

Inequalities (|3.2|),(|3l^) greatly simplify the problem. Firstly, ( |373| ) implies that the decay 
takes place when the Wannier pair is far outside the atomic radius. As a result the 
electrons e^e/j have no influence on the process of ejecting the electron e^. Secondly, 
inequality (|3.2|) guarantees that the interaction between the Wannier pair and the 
autoionized electron is small, and we can neglect it. As a result the matrix element 
(^f|VA|^j) reduces to 

(^fl^Al^j) ^ ^ft\{^2\t)- (3-4) 
Here Ff is the partial width of the autoionizing state, which makes \/Tf proportional 
to the matrix element responsible for the autoionization decay. It is determined by the 
processes which take place strictly inside the atomic particle and, as mentioned, are 
not influenced by PCI. The effects of PCI are taken care of in the overlapping integrals 
(V'2|V'j) between the wave functions which describe the states of the Wannier pair. We 
denote by {ip2\ the flnal-state wave- function for the Wannier pair and by {ipj) the full 
set of the intermediate-state wave functions for the pair. Substituting (|3.4|) in (|3.1|) we 
can present the matrix element M in the following simple form 

M = Jl\{^2\^Pi). (3.5) 
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Here the wave function {ipi) is defined as 

\^,) = {^{A^*)\GiX\r/,\V.oi\^d, (3.6) 

where (\[^(A+*)| is the wave function of the resonant state of the ion. It projects the 
many-electron wave function Koil^i) onto the states of the Wannier electrons. The 
function depends on the coordinates of the two electrons. Summation over the full 
set of intermediate states for the pair is included in the two-particle Green function 

M+) ( I ^j(ri,r2)^;(r;,r^2) 

Equation ( |3.5|) states that the influence of PCI on the Wannier pair can be described 
with the help of the overlapping integral. This result is very close in spirit to the shake- 
down model |jl3|, |I4| which originally was developed for the case when one near-threshold 
electron is exposed to the variation of the charge of the atomic particle caused by the 
autoionization. In the case under consideration there are two near-threshold electrons, 
but the idea that overlapping of the wave functions gives the correct description of PCI 
remains the same. 

The wave functions in the matrix element (p.5|) satisfy the Schroedinger equations 



£ + z^-i7i)|^i) = |g), (3.8) 
{E-H2)\^2) =0. (3.9) 

Here Hi and H2 are the Hamiltonians which describe propagation of the Wannier pair 
in the intermediate and final states 

Hn = -^(A„ + A^) + - Z„ i^ + Pj. (3.10) 

An index n = 1,2 here labels the intermediate and the final states in which the 
Coulomb charge of the atomic particle is different. Generally speaking there are also 
the short-range potentials which account for the potential of the atomic particle in the 
intermediate or final states. However, the influence of the short-range potential on PCI 
is small, and these potentials are neglected in ( p.lO| ) and the following consideration. 
The energy e + iV /2 of the pair in the state iV'i) is complex. Its real part e is equal to 
the above-threshold energy, its imaginary part V /2 appears due to autoionization decay. 
The function \il)2) describes the Wannier pair in the final state with real energy E which 
satisfies the energy conservation law E^oi + E{A) = E + e^ + E{A^^), where E{K) is the 
energy of the atom, E{A^^) is the energy of the ion, is the energy of the autoionized 
electron, and E^ox is the collision energy. 

According to ( |3.(j| ),( P^ there appears a source in the right-hand side of the 
Schroedinger equation (3.8) 

|g) = (v[/(A+*)|Koi|^.) (3.11) 
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which is a function of the coordinates of the pair. For large separations from the ion 
this function becomes proportional to the amplitude of inelastic scattering and inversely 
proportional to separation. It is important that the effects of PCI take place far outside 
the atomic particle 0, where the source is small (rQ,r/3|(5) ^ 0, rQ,r/3 ^ 1. Therefore 
it is sufficient for our purposes to solve only the homogeneous Schroedinger equation for 
the function \ipi) 

[e + i^-H^^\^^) = (3.12) 

in the region far outside the atom and then normalize this function to make it 
proportional to the amplitude of inelastic scattering, thus taking into account the 
influence of the source. 

The given consideration shows that PCI for the Wannier pair in the vicinity of the 
threshold is described by the overlapping integral of the wave functions describing the 
pair in the intermediate and final states ( |3.5| ). Notice that this possibility follows from 
conditions ( ^.2[ ),( P^ ). In this sense it is rather a theoretical result then a model. 



4. PCI for the Wannier pair 

Let us calculate the overlapping integral (p.5|) using the wave functions for the Wannier 
problem evaluated in Section |^. The model suggested in Section ^ describes the state 
of the Wannier pair by its energy and the quantum number x- The amplitude of the 
inelastic collision depends on these quantum numbers 

^(x) = (e,x;^(A+*)|Koi|^i). (4.1) 

The matrix element in the right-hand side describes the events which take place in the 
vicinity of the atom. The small width F of the autoionizing state does not manifest 
itself for these small separations. That is why we can neglect it, supposing that the 
wave function {e, x\ in ( |4.1|) describes a real physical state of the pair with the real 
energy e. 

The solution of (|3.12|) can be presented in the following form 



l^i) = J [i^k,x'{r,r])]^Aix')dx'- (4-2) 
Here the wave function ipk,x' ^7) defined in ( |2.62|) describes a state of the Wannier pair 



['^k,x'{r,v)]i = \j B{x) exp[-iS(r, rj; k, xM^t^\r] w{x')) 



(4.3) 



The square brackets with the index 1 are used to remind that the Wannier pair is 
in the intermediate state in which the Coulomb charge of the ion is Z = Zi = 1. 
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Asymptotic condition ( |2.71| ) guarantees that the coefficients A{x') i^i identical 
to the amphtude of inelastic scattering. 

It is important that dependence of the amplitude A{x') oii x' can be established 
explicitly. Remember that ( |2.62| ) expresses the properties of the Wannier problem in 
terms of the effective single-particle problem. This allows one to present the amplitude 
of the collision as the following product 

A{x') = exp[z6{x')]Ao{x') ^ exp[z5(x')]A(x)- (4-4) 



Here S{x') is the phase shift which was found in (^77^ ,(^77^ ) ,(^7R]. The quantity 



^o(x') is a matrix element of inelastic collision. Representation ( [4.4|) has a form which 
is conventional in the processes with single-particle final states [ill] . The matrix element 
Ao{x') originates from integration in the vicinity of the atom. In contrast, the phase 
shift 6{x) arises due to the events which take place far outside the atom in the Coulomb 
zone. Thus the two quantities Ao{x') and S{x') originate from quite different distances. 
This makes them depend differently on x'- The phase shift exhibits strong variation 
Hx') ~ x'^/^! while the amplitude Ao{x') varies smoothly ^o(x') ~ x'^- order to 
verify this statement one should remember that in the semiclassical region all quantities 
depend on the classical action ~ / dr^ k? + w(x')/r. For large r, r ^ fc^, the action 
has a term ~ x''^ /k which strongly varies with x', while for small r, r -C fc^, there is 
only a smooth behaviour ~ x'^ <^ 1. More general verification of this property can 
be fulffiled using the Landau-Smorodinsky approach to the problem of a low-energy 
particle influenced by a long-range Coulomb field as well as a short-range potential. It 
is discussed in detail in 



The smooth behaviour of the matrix element Ao{x') with x' justifies the last equality 
in Using it we simplify 



= Aoix) j [exp[z5(x')]^fc,x'(^,^)]i^^x'. (4.5) 

The momentum k in the state is a complex number due to the finite width of the 
autoionizing state 



k = V2e ir ~ V2e + i — (4.6) 



Inequality ( p.3|) permits one to make expansion in the last equality in ( [4.6| ). It is 



important to keep the width in this expansion because the wave function ( [4.2| ) will be 
used to calculate the effects of PCI which take place for large separations. The positive 
sign of the imaginary part of k ensures that the wave function decreases exponentially 
['4^k,x{''^yV)]i ~ exp(—Im /cr) with separation. Thus the events which result in PCI are 
localized in the finite region. Similar localization of the wave function is well known 
M for the case when one electron is influenced by PCI. We see that localization in the 
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vicinity of the atom 

/ 1 2V2i 

r<^^^^n (4.7) 

remains true for the Wannier pair as well. The radius of localization ri has a clear 
physical meaning being proportional to classical distances ri(t) ~ r2(t) ^ ri/(2\/2) 
which separate Wannier electrons from the atom at the moment t = l/F of the decay. 

The wave function {iIj2\ in the overlapping integral ( |3.5D describes the real final 
state of the Wannier pair with the real energy E. Generally speaking this state 
should be described by a linear combination of outgoing and incoming waves. For our 
purposes it is necessary to take into account the convergent wave because it gives the 
major contribution to the overlapping integral and therefore correctly describes PCI. 
Really, the convergent wave function behaves like ~ exp{—ipr), p = \/2e therefore 
the overlapping integral with the divergent wave function ~ exp(i/cr) can be roughly 
estimated as ~ Xjik — p) ~ l/F describing the resonance. In contrast, the divergent 
wave does not contribute to the resonance. This shows that the final-state wave function 
|^/'2) can be taken in the form (|2.62| ) 



1^2) = \^vA°^^^\ 



\jB{x) exp[zS(r, 77; p, x)]</>i, Hr;w(x)) 



(4. 



The square brackets with the index 2 are to remind that all the quantities in the square 
brackets should be calculated for the final state in which the Coulomb charge of the 
ion is Z = Z2 = 2 and the momentum of the pair is p = ^1E. Notation introduced in 
,([4. 81) will be used below as well. All quantities labeled by an index 1 are considered 



for the intermediate state of the Wannier pair. In this state the ion charge is Z = Zi = 1, 
the energy of the pair is e+iF/2, the corresponding momentum is k ([4.6|) . The quantities 
labeled by an index 2 describe the final state with Z = Z2 = 2, the final-state energy of 
the pair E and the corresponding momentum p = \JlE. 

The explicit form of the wave functions ( [4.5| , (|4.8| ) permits to present the matrix 
element describing PCI ( p.5|) in the following final form 



M = ^Ff A(x)52i, (4.9) 

where 

^21 = / [V^,,,(r,r/)]2[V^.,,.(r,r^)]iexp[^5i(x')]^^!^- (4.10) 
J Zn 

The above given consideration shows that this representation has a simple physical 

meaning. The fast autoionization decay results in an abrupt transition of the Wannier 

pair from the intermediate state to the final state. This transition is described by the 

overlapping integral of the wave functions for the intermediate and final states 

/d/Tidv 
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The integration over dx' in ( [4 .101 ) reflects the fact that the intermediate state includes 



the coherent superposition of the states of the pair with different x' excited due to 
inelastic collision. The factor exp[i5i(x')] describes fast variation of the amplitude of 
the excitation with x'. 

Integration in ( [4.10|) over drj is Gaussian, therefore it can be fulfilled analytically, 
allowing the presentation of the amplitude S21 as a two-dimensional integral which can 
be easily handled numerically. 

Knowing the matrix element M one can present the cross-section of the reaction 
(PD,(pri) in the form 

dcy^ci{E) = ^R^,,{E,e)da^{e). (4.12) 

The cross-section da^ciiE) describes the full reaction in which PCI plays a role. The 
cross-section dao^e) describes only reaction ( p..l| ) in which the role of PCI is neglected. 
The factor 



r 



i?pei(^,£) = --|52l|' (4.13) 

/vr 

presents the influence of PCI in the explicit form described by the amplitude S21 in 
( [4.10| ). This factor has a resonant nature sharply depending both on the above-threshold 



energy e and the final-state energy of the pair E. 

Equations ( 4.10 ),( 4.12 ),( 4T3|) explicitly describe the influence of PCI on the cross- 



section. The physical consequences of these formulae are discussed in the next Section 
i 

5. Manifestations of PCI 

5.1. Lorentz line 

Consider first the trivial situation when one neglects the effects of PCI supposing that 
the Coulomb charge of the ion in the intermediate and final state remains the same 
Zi = Z^. Then the wave functions marked by symbols 1 and 2 in (|4.2| ),( ^^ differ only by 
the momenta k and p. Notice that the coefficients A[x\ B{x), C{x) in the expansion of 
S(r, rj] p, x) (|2.571 ) are smooth functions. Therefore using inequality (|3.3| ) we can assume 
that they do not vary much when the energy parameter e in ( p.l6|) takes either the value 
e = E, OT e = e, i.e. A{x)\,=e ^ A{x)\,=e, B{x)\,=e ^ B{x)\,=e, C{x)\,=e ^ C{x)\^=e. 
The last equality greatly simplifies integration over drj in ( [4.10|) 

drj 1 



J exp[iB{x){x - X')v]^ = '^(^^^{x - x')- 



Here 5d{x ~ x') is the Dirac delta-function, which is eliminated by the following 
integration over dx' in ( [4.1UD . As a result the amplitude can be presented as the simple 
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one-dimensional integral 

S21 = exp[i5(x)] / (\t^{r]w{x)Wk \^'^w{x))^'^ (5-1) 





with the Coulomb wave functions in the field of the same effective Coulomb charge 
Evaluating it with the help of (|2.64|) and substituting into ( [4.13|) one finds that 
the resonant factor RpciiE,e) reduces to the Lorentz line Rii{E,e) 

T/2TT 

{E-ey + ryr 



Rpci{E,e) Rn{E,e) = J,, , (5.2) 



as it should be when PCI is neglected. Notice that this result justifies the normalization 
of the wave functions in ( p.64|) . 



5.2. Eikonal region 

Let us examine now the infiuence of PCI on the Wannier pair. The integration over 
dr in (^4.101) is localized in the region r < n ( [4.7|) . It is essential how the localization 



radius ri is related to the Coulomb radius Vc (|2.34]) . Let us examine the situation when 
condition 



is fulfilled. It is valid for sufficiently high energies E, e. In this case the main 
contribution to the integral ( [4. 101) comes from the region where asymptotic relations 
(2.68), (2.69), (2.70) are fulfilled. They imply that Bi{r,k) ~ B2{r,p) ~ 1 and 
Ci(r, k) ~ C2{r,p) ~ 0. These equalities show that the integration over drjdx' in (|4.10| ) 
gives the trivial result 

exp[i{x - x')^]^^^^^ = 1- (5-4) 

Using it we find that the matrix element 5*21 is simplified to be proportional to a one- 
dimensional overlapping integral 

5'2i = exp[z(5i(x)](</)2|(/)i), (5.5) 

(02|0i)= J^^ 4-\r;w2ix) W\r;wiix))dr. (5.6) 

Both wave functions here describe the effective single-particle Coulomb problem. This 
problem in the intermediate state is governed by the effective charge Wi{x) defined by 
( p.6| ),(2.7),(2.8) with the ion charge Z = Zi = 1. In the final state the effective charge 
W2{x) depends on the final-state ion charge Z = Z2 = 2. 

From (5.5),( ^?TB| ) we find the following simple representation for the resonant factor 

Rp,i{E,e) ~ R,{E,e) = ^|(02|0i)|'. (5.7) 
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The factor Rc{E,e) defined in (5.6) depends on the overlapping integral between the 
wave functions describing the effective single-particle Coulomb problem. It is remarkable 
that this integral proves to be similar to the quantity which was well-known previously 
in the traditional PCI problem when a single low-energy electron is influenced by the 



products of the autoionization decay. In this case the shake-down model |TB|, [T^, ^ 
gives correct description of PCI. The shake-down model is based on an overlapping 
integral with single-electron wave functions, whose structure is very close to the integral 
in (5.6), (^. 

We come to a very interesting conclusion. Our knowledge of PCI properties of 
the traditional situation when there is only one near-threshold electron enables us to 
predict manifestations of PCI when the Wannier pair is involved. Conveniently we do 
not even need to fulfill explicitly calculation of the overlapping integral (5.5) which 
appears in { ^77\j . Inequahty (|5.3| ) guarantees that after the decay the Wannier pair 
is located so far away from the atom that neither the Coulomb charge of the ion nor 
repulsion between the electrons can significantly change the classical trajectory. Thus 
the pair moves almost without acceleration, which greatly simplify behaviour of the 
wave functions describing this motion. A similar case was well studied in the single- 
electron problem. It was shown that the eikonal approximation and the shake-down 
model give the same results when conditions of applicability for the eikonal theory are 
fulfilled. Borrowing this result from the true single-electron problem and applying it to 
our effective single-particle problem we conclude that the eikonal theory proves to be 
correct for the Wannier pair if condition ( p. 3] ) is fulfilled. Notice that condition ( |5.3| ) 
can be rewritten in the standard in the eikonal approximation form e^/^ ^ P. Of course 
one could anticipate from the very beginning that when this usual condition is valid, 
then the eikonal theory should give the correct description of PCI for the Wannier pair. 
However, the peculiarity of the Wannier problem certainly appeals for a steady basis for 
such a statement. That is why it is important that we have been able to demonstrate 
validity of the eikonal theory by direct calculations. 

Recognizing that the eikonal approximation is valid, we conclude that PCI for the 
Wannier pair should have the following manifestations. 

1. The line in the autoionization spectrum is shifted by PCI. The shift is equal to |15, ^ 

AE= - (5.8) 

e = l + (5.9) 

V2i V, v^ 

Here Vq, = v^ = ^/e. The right-hand side in (5.9) is presented in the form which is 
common in the eikonal approximation. 

2. The line changes the shape, becoming broader, lower, and acquiring the "left 
shoulder" which makes it asymmetrical. All these variations are the more significant 
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the larger is the basic parameter ^. These properties of the hne are described by the 
profile ^ I 

R^,i{E,e) ^ Reii.{E,s) = ^^^-^^^L__k,^^^E-e)/T,0, (5.10) 
keikHE — e) /T , ^) = , ^ exp [— 2^ arctan- 



sinhvre ^ r/2 

3. In spite of strong variations in the shape and position the intensity of the line remains 
constant, it is not influenced by PCI, 

I Rp,iiE,e)dE = l. (5.11) 

This property is essential for the Wannier problem where intensity is one of the main 
characteristics of the process. Taking into account that the cross-section of the reaction 
( |1.1|) satisfies the power law (|1.2| ) with the index fi = fii = 1.127 we find from 



|1.2|) ( [4.12D , (|5.11|) that this law remains valid for the total intensity of the line for reaction 

dE ao e^"' (5.12) 



(0),(e: 

d(Tpci{E) 



a 



dE 

considered as a function of the above-threshold energy. 

Up to this point our consideration was restricted by two particular assumptions. 



Firstly, we considered the case of large velocities of the autoionized electron ( |3.2|) which 
allows one to assume that it does not play a role in PCI. Secondly, we discussed the 
opposite directions for the Wannier electrons neglecting fluctuations in their angular 
distribution. Once we have established that the eikonal theory is correct we can remove 
both these simplifications and include into our scheme the case of sufficiently low-energy 
autoionization and arbitrary directions for the Wannier electrons. As usual in the eikonal 
approximation it is sufficient to redefine the parameter ^ which has the following more 
general form 

,.L.L-^-±. (5.13) 

Here Vq,^ = Vq, — v^, v^^ = — v^. With this parameter ^ (5. 8), (5. 10) describe the 
line which depends on the absolute values as well as directions of the velocities of all 
three electrons. 

5.3. Coulomb region 



Consider the case when the localization radius ri (^.7|) is smaller than the Coulomb 



radius (|2.34|) 



n = -jr- < = (5.14) 



25 



This inequality holds for sufficiently low energies E, e. This inequality means that the 
decay takes place when the pair is so close to the ion that the Coulomb potential has a 
strong influence on trajectories of the Wannier pair. 

Let us examine the behaviour of the integrand in the amplitude 5*21 in ( [4 .101 ) as 
a function of x'- Remember that in the region r <C Tc the asymptotic relations 
(2. 65), (2. 66), and (2.67) for the coefficients governing the action S are valid. Notice 
ffist of all that the quadratic in x' terms in the exponent arising from the phase 
shift A6i{x') cind from the term [T,{r,ri;p,x')]i cancel each other due to (2.65),( P?7^ ). 
Another simplification comes from the fact that for small distances the Coulomb function 
smoothly depends on x' exhibiting behaviour ~ x'^, see discussion after ( [4.4|) . Therefore 
one can suppose that 



[exp[z5c(x')]0fc'^^(^;^(x'))]i - [exp[i5c(x)]0r^('^;«^(x))]i- 
These properties of the integrand allow analytical integration over dx' in ( |4.1(j| ) 

dx' 1 



f exp{-i[B{x)]ix'v)^ 



(5.15) 



The delta-function is eliminated by the subsequent integration over di]. As a result one 
finds the following simple representation for the amplitude 

= exp[z4,i(x)](02|(52/5i)i/2|0i). (5.16) 

Here the matrix element is defined as 



t\r; W2{xmB{^)]2/[B{x)]^Y/'4-\r; W^{x))dr. 



Using ( |5.16| ) we find the following representation for the resonant factor ( [4. 131 ) 



r 

2^ 



B2/B^\(P,) 



(5.17) 



(5.1^ 



This expression can be simplified even further. The functions [-B(a;)]i and [5 (a;)] 2 
smoothly depend on the coordinate r. In contrast, the Coulomb wave functions |02) 
oscillate. Calculating the integral with sufficiently fast oscillating functions one can use 
the saddle-point method. The semiclassical phases of the wave functions \4>i), \4>2) are 
/ ^2{e + iT/2 + iyi(x)/r)(ir and / \J2{E + W2{x) l''^)dT respectively. Their difference 

dr. (5.19) 



IS 



<l>(r) 

The saddle point r 

E + 



sp 



W2{X) 



sp 



2{e + iT/2 + Wi{x)lr) 



is to satisfy equation $'(r, 

^ '"sp 



spy 



2{E + W2{x)lr) 
= yielding 



(5.20) 
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Solving this equation one finds 



' sp 



(5.21) 

is large. This fact 



e-E + iV/2 

The energies E,e are low (|5.14| ), therefore the radius r = r^^ 
justifies both the applicability of the semiclassical estimation for the phases in (|5.19| ) 
and the validity of the saddle-point approximation. These arguments demonstrate that 
calculating the matrix element ( |5.17 ) we can suppose that the main contribution to the 
integral comes from the vicinity of the saddle point (|5.21| ). This allows one to evaluate 
smooth functions [i?(x)]2 at this point and take them out of integration. As a 

result the resonant factor given by ( |5.18 ) can be simplified to the following final form 
which was first announced in ( |1.5| ) as the major result of this paper 

(5.22) 



R^,,{E,e) = K^{E,e)R,{E,e) 



where the factor 



(5.23) 



depends on the overlapping integral (02|0i). Notice that this integral can be presented 
in an analytical form if the integration is pushed one step further using the saddle-point 
method, but for our purposes it is sufficient to keep it in general form. 

We met the factor (|5.23|) before when discussed the eikonal approximation, see (|5.7|) . 
It was argued there that the factor Rc{E, e) describes those properties of PCI which were 
previously well known from traditional studies of PCI with one low-energy electron. In 
contrast, the factor 

K„{E e) = ^^^^^'^P-* 

which appears in (|5.22| ) is a new quantity specific for PCI with Wannier electrons, having 
no analogue in PCI with one low-energy electron. Due to this reason it is called the 
Wannier factor. Remember that the function B{x) defined by (2.51),(2.26), (|2.16|) gives 
normalization for the wave function of Wannier electrons (|2.62| ). The found Wannier 
factor depends on a ratio of these normalization functions for the intermediate and final 
states. The coordinates of normalization functions, in accordance with ( p.l6| ) are defined 



(5.24) 
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X2. 



sp 



sp 



1 + 



1/2 



+ 1 



er. 



sp 



1/2 



Er. 



sp 



(5.25) 



(5.26) 



Traditionally the interest in the Wannier problem has been inspired by the near- 
threshold power-law of the cross-section ( [1.21 ). It originates from the asymptotic 
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behaviour of the normahzation coefficient B{x) for small separations when Er <^ 1, see 
(2.66). Thus traditional studies of the Wannier problem can be considered as a probing 
of the normalization coefficient B{x) in the particular asymptotic region well inside 
the Coulomb zone. PCI opens an interesting opportunity to investigate this coefficient 
in a much wider area, not necessarily restricted by the Coulomb zone. According to 
(5. 25), (5. 26) variation of the above-threshold energy e and the final-state energy E 
leads to variation of Xi^gp, X2,sp in a broad area resulting in significant variation of 
normalization coefficients [S(a;)]i and [i?(x)]2. 

The result obtained ( |5.22D permits one to address the question of the threshold 
behaviour of the cross-section which is of particular interest for the Wannier problem. 
When the energy of the Wannier pair is low E Q then one finds from (|5.21|) ,(5.26) 
that ~ iW^ix) ~ W^i(x))/^ aiid X2,sp ~ ^I^J~Ei oo. The later condition allows the 
use of asymptotic relation (2.66) which shows that B2[x2,sp) ~ E^'^ ^ ^2 = ^^2/2. From 
( ^.24| ) one finds the same estimate for the Wannier factor 

R^,i{E,e)^E>^\ (5.27) 

Equation ( |5.22| ) gives the same estimate for the resonant factor Rpci{E,e) ~ E'^^ 
resulting in the threshold law for the differential cross-section as a function of the energy 
of the pair E 

( da\ 



dE 

/ 1 



E-^O. (5.28) 

pa 

Similarly one can estimate the behaviour of the resonant factor on the above-threshold 
energy Rpci{E,e) ~ £:^^\ = Vxjl. Remembering that the cross-section of the 
exhibits the usual Wannier power-type behaviour o"o ~ £ — > as a function 

of the above-threshold energy we find that the cross-section for the combined process 
( prTD , (|1.4|) does not depend on the above-threshold energy e 

f da\ 



dE / 



~ const, e ^0. (5.29) 

pet 

We see that PCI has a dramatic effect on the threshold behaviour of the cross-section. 
Firstly, the PCI results in the variation of the exponent. The exponent in ( 5.28| ) 
corresponds to the final-state charge of the ion Z = Z2 = 2 which gives fi2 = 1-056. If the 
Auger decay is impossible, then the cross-section of the reaction ( p. . 1|) is described by the 
usual Wannier power law ( |1.2| ) which corresponds to the charge of the ion Z = Zi = 1 
resulting in /ii = 1.127. Secondly, the cross-section does not depend on the above- 
threshold energy e. These results agree qualitatively with the simple physical picture. 
If above threshold energy is low then the decay takes place when the pair is so close 
to the ion that its separation can be neglected compared with the large radius of the 
Coulomb zone. This means that the propagation of the pair in the intermediate state 
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is insignificant, which makes insignificant as well the very existence of this intermediate 
state. That is why the final answer for the Wannier exponent should not depend on 
parameters governing the intermediate state. The fact that ( |5.22| ) reproduces this result 
can be considered as a qualitative verification of this formula. 

Generally speaking, one could contemplate a possibility to measure the power-type 
behaviour ( p.28| ) . To this end one should fix the above-threshold energy and measure the 
"left shoulder" of the resonance profile in the near-threshold region versus the energy 
of the Wannier pair. However this project meets a difficulty. It can be verified that 
in order to distinguish the exponent fi2 from fii in the power-type behaviour of the 
cross-section (|5.28|) the left-hand side in inequality (|5.14| ) should be really small, say, 
less than < 10^^ putting a severe restriction on E,e. 

The most important result of this Subsection is ( ^.22| ) which describes the infiuence 
of PCI on the resonance profile. This equation was evaluated assuming that energies E, e 
are low enough to satisfy condition ( ^.14|) . Notice, however, that the equation remains 
correct for much higher energies satisfying the condition of applicability of the eikonal 
approximation (|5.3|) which is opposite to the low-energy limit. To see this consider 
the eikonal region (|5.3| ) in which — e| ^ F^/^. For these energies (|5.21| ) 

results in Tgp ~ F ^^"^ and therefore Evgp ~ £T^sp ^ 1- Equations (5. 25), (5. 26) show 
that Xi^sp ~ X2^sp — 1 permitting one to use asymptotic condition (2.69) which shows 
that the Wannier factor is trivial K„{E,e) ^ 1. Therefore, in the eikonal region (|5.22| ) 
yields Rpa{E,e) ~ Rc{E,e) which according to ( p.7| ) is the correct result. The fact that 
( p.22| ) proves to be correct in both the eikonal region r\ <^ and the Coulomb region 
r\ ^ Tc may be considered as an indication that it should give reasonable results in the 
intermediate region r\ ~ as well. 



6. Electron impact in the vicinity of Ne K-shell 

In order to illustrate the validity of different approximations developed above for 
description of PCI let us consider the example studied experimentally in in which 
the electron impact ionization of the Ne atom in the vicinity of the K-shell 

e + Ne ^ e + e + Ne+(ls"i) (6.1) 

is followed by the KLL Auger decay 

Ne+(ls-^) ^ e + Ne2+(2p^2)_ (g 2) 

The distortion of resonance profiles by PCI is illustrated in figure 1 for above- 
threshold energies e = 2,3,5 and 10 eV. The width of the K-vacancy is F = 270 meV 
p!6| . The profiles are calculated for the symmetric geometry of Wannier electrons when 



they move in opposite directions possessing equal energies. Notice that experimentally 
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|T|] the spectra were measured versus the energy of the Auger electron when neither 
energy nor angular distributions of Wannier electrons were fixed. 

The full curve in figure 1 presents the results of calculations based on the most 
advanced equation ( p. 221 ) which, as was argued above is valid for both near-threshold 
and high energy regions. The factor Rc{E,e) = r/(27r)|((^2|(^i)P is calculated by direct 
numerical integration of the overlapping integral with the Coulomb wave functions 
describing the effective single-particle problem. Remember that this problem gives the 
simplest description of the Wannier pair. Normalization of the Coulomb wave functions 
is given in (|2.64|) . For numerical calculations it is more convenient to work with the real 
Coulomb wave function 02 (^) taking a linear combination of convergent and divergent 
waves, rather then only the convergent wave. This is possible because the divergent 
wave in 02 (^) does not contribute to the resonance, see discussion before (|4.8| ). The 
wave function 0i(r) has to be taken as a divergent wave localized in the vicinity of the 
atom ip. The Wannier factor K„{E,e) is found from (2. 51), (2. 26). 

Full dotes in figure 1 present calculation based on ( |5.18|) . The perfect agreement 
between the full lines and the full circles shows the high accuracy of the approximation 
which was used in order to present the resonant factor in the final form ( ^.22|) as a 
product of two different factors. 

Dotted lines in figure 1 present the simple eikonal approximation (5.10) which should 
be correct only for high energies. Results presented show that below 6 = 5 eV the eikonal 
approximation fails, while for high above-threshold energies e > 10 eV it produces quite 
reliable results. Notice that experimental results of belong to the region e > 10 eV 
where we have just verified the eikonal approximation to be valid. 

The dotted line in figure 1(a) shows the Coulomb factor Rc{E, e) which is known to 
produce good results for PCI with a single low-energy electron. Figure (la) shows that 
for the considered situation with two low-energy electrons it proves to be incorrect for low 
energies. This demonstrates that the Wannier factor K„{E,e) plays a very important 
role. The Wannier factor exhibits quite an interesting variation as demonstrated in figure 
2 where it is presented versus the energy of the pair E for the fixed above-threshold 
energy e = 2 eV. In order to reveal the asymptotic low-energy behaviour ~ E'^'^ , E 
predicted by (|5.28|) the Wannier factor is shown as K^{E,e)/E. 

The PCI shift found from calculations of the PCI profiles is shown in figure 3 for the 
symmetrical configuration of Wannier electrons. The full curve gives the results obtained 
from ( 5.22|) . The dotted line shows the prediction of the eikonal approximation. The full 
dottes in figure 3 present the results extracted from the profiles calculated in 0. These 
results effectively include different possible configurations of Wannier electrons in the 
final state. One can expect that the symmetrical configuration of the Wannier electrons 
gives large contribution to the total probability. Therefore the results calculated for 
the symmetrical configuration should reasonably agree with calculations in which all 
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continuum states of Wannier electrons are included. Notice that the theory developed 
above and the approach of are very different in technique used. Nevertheless there 
is a clear resemblance in several important basic features. Firstly, in both works PCI is 
described with the help of overlapping integrals between wave functions describing the 
Wannier pair. Secondly, both works rely upon the semiclassical approximation. Having 
in mind these basic similarities one could expect that the two approaches should give 
similar results. This hope proves be correct for high energies, see figure l{d) for 10 
eV above the threshold. There is, however, some discrepancy below 10 eV. It can be 
verified comparing figure 1 with the results presented in figure of |^ that distinctions 
exist in profile shapes as well. 



7. Discussion of the main results 

This paper develops the following ideas and results. 

Lit is emphasized that the propagation of two Wannier electrons can be described 
by an effective single-particle Coulomb problem with corrections which take into account 
more subtle three-body properties. Correspondingly the wave function of the Wannier 



pair ( p.62|) is presented as the product of the Coulomb wave function describing the 
effective single-particle Coulomb problem, the additional phase factor exp(zE), and the 
normalization coefficient B (2.51) which plays a very important role in our consideration. 

2. The complex problem ( |1 . 1| ) , ( p^ ) with three electrons in the final state is reduced 
to a more simple problem in which only two Wannier electrons play a role. The key 
point is an assumption that the autoionization energy is so high that the autoionized 
electron cannot interact strongly with the other two electrons. The dynamical properties 
responsible for PCI are shown to be described by the overlapping integral (|4.10|) between 
the wave functions of Wannier electrons in the final and intermediate states. 



3. The resonant factor is found in a very simple form (|5.22|) which depends on 



two factors. One of them Rc describes the overlapping integral between the Coulomb 
wave functions for the effective single-particle problem. A similar integral is known to 
describe PCI when there is only one near-threshold electron. More specific properties 
of the considered problem are described by the Wannier factor The factorization 
of the resonant factor into a product of and Rc has a clear physical reason. There 
are two scales in the problem. The effective single-particle Coulomb problem describes 
strong variations of the wave function with a typical radius Vgi^gic ~ In contrast, 
the quantities S, B describing the specific properties of the three-body problem vary 
with a typical distance equal to the radius of the Coulomb zone r3_body ~ ^oi low 

energies these radii are different rsmgie -C r3_body It is well known that when there are 
two different scales in a problem then one has to expect an amplitude and a probability 
to be presented as products of quantities describing events which happen in different 
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scales, as it happens in our case. 

Equation ( p. 22 ) was verified for low energies, i.e. in the Coulomb region, and for 
high energies in the eikonal zone. For intermediate region one can hope that it should 
give reasonable results. (However, validity of the later assumption should be clarified 
by numerical integration of the matrix element ( [4.10[ ).) 

The Wannier factor is shown to depend on the normalization coefficient B of the 
wave function of the Wannier pair, see ( p. 241 ). The distance at which this coefficient is 
to be taken depends on the above-threshold energy as well as on the energy of the pair 
in the final state, see ( ^.21| ). It is remarkable that the same coefficient considered for 



small separations (2.66) governs the Wannier power law ( |1.2|) , and therefore it has been 
under thorough experimental investigation for a long period of time. Equation (|5.22| ) 
shows that there is a new possibility to experimentally investigate the coefficient B in 
a wide range of distances. 

4. It is demonstrated that for low energies of the pair E ^ the Wannier power 
law has the form (|5.27|) with the exponent governed by the ion charge in the final state. 
In contrast for higher energies in the eikonal region the power law has a form of (|5.12| ) 
depending on the above-threshold energy e with the exponent governed by the ion charge 
in the intermediate state. 

5. Our consideration was restricted by the one-dimensional model formulated in 
Section ^ However, one can hope that the approach developed is more general than 
the model itself. There are clear physical reasons for this. Firstly, PCI as well as 
the Wannier problem manifest themselves for large separations where the Coulomb 
potential play an important role, while the more rapidly decreasing orbital barrier does 
not. Therefore neglect of this later potential seems reasonable. Secondly, the arguments 
considered above show that the factorization of the resonant factor into the product of 
the Coulomb factor and the Wannier factor is a general property caused by existence 
of the two different scales in the problem. Therefore it should remain valid when the 
angular variables are included as well. 
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Figure captions 



Figure 1. Electron impact with excitation of K-shell in Ne. Calculated profiles of 
Auger lines distorted by PCI R{E, e) versus the energy of the Wannier pair E for a 
fixed above-threshold energy e. The width of the K- vacancy is 270 nieV Figures 
(a),(&),(c),and (d) present results for e = 2, 3, 5 and 10 eV. The symmetric configuration 
of Wannier electrons is considered, they move in opposite directions possessing equal 

energies. — • — the non-distorted Lorentz line, prediction of the most advanced 

equation ( ^.22 ), the eikonal approximation Rcik{E,e) (5.10), • prediction 

of equation ( 5.18 ), the Coulomb factor Rc{E,e) ( 5.23| ). The perfect agreement 

between the full line and full-circled one shows the high accuracy of the approximation 
which permits to present the resonant factor in the final form ( 5.22 ). Figure (d) shows 
that the eikonal approximation is valid above 10 eV. The strong deviation of the dashed 
line from the full one in figure (a) indicates an important role played by the Wannier 
factor Kv,{E,e) for low energies. 



Figure 2. The same reaction as in figure 1. the ratio of the Wannier factor 

K^{E, e) to the energy of the pair E versus the energy E for the fixed above-threshold 
energy e — 2 eV, low-energy asymptotic behaviour ( ^.2^ ). 



Figure 3. The same reaction as in figure 1. The calculated PCI shift versus the 

above-threshold energy: prediction based on equation ( 5.22| ), the eikonal 

approximation (5.10), ©results extracted from profiles calculated in ||]. 
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